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Abstract 

The problem of escape of a Brownian particle in a cusp-shaped metastable potential is of special 
importance in nonadiabatic and weakly-adiabatic rate theory for electron transfer (ET) reactions. 
Especially for the weakly-adiabatic reactions, the reaction follows an adiabaticity criterion in the 
presence of a sharp barrier. In contrast to the non-adiabatic case, the ET kinetics can be, however 
considerably influenced by the medium dynamics. 

In this paper, the problem of the escape time over a dichotomously fluctuating cusp barrier is 
discussed with its relevance to the high temperature ET reactions in condensed media. 
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I. INTRODUCTION 



Mechanism of the electron transfer (ET) in condensed and biological media goes beyond 
universal nonadiabatic approach of the Marcus theory. In particular, relax- 

ation properties of medium may slow down the overall ET kinetics and lead to an adiabatic 
dynamics. (?] An excess electron appearing in the medium introduces local fluctuations of 
polarization, that in turn contribute to the change of Gibbs energy. Equilibration of those 
fluctuations leads to a new state with a localized position of a charge. In chemical reactions, 
the electron may change its location passing from a donoring to an accepting molecule, giv- 
ing rise to the same scenario of Gibbs energy changes that allows to discriminate between the 
(equilibrium) states "before" and "after" the transfer (see Fig. [IJ. The free energy surfaces 
for "reactants" and "products" are usually multidimensional functions which intersect at the 
transition point. The deviation from it, or the Gibbs energy change, can be calculated from 
the reversible work done along the path that forms that state, so that by use of a simple 
thermodynamic argument, one is able to associate a change in the Gibbs energy with the 
change of multicomponent "reaction coordinate" that describes a response of the system to 
the instantaneous transfer of a charge from one site to another. 

ET reactions involve both classical and quantum degrees of freedom. Quantum effects 
are mostly related to electronic degrees of freedom. Because of the mass difference between 
electrons and nuclei, it is frequently assumed that the electrons follow the nuclear motion 
adiabatically (Born-Oppenheimer approximation). The interaction between two different 
electronic states results in a splitting energy A. The reaction from reactants to products 
is then mediated by an interplay of two parameters: time of charge fluctuations between 
the two neighbouring electronic states and a typical time within which the nuclear reaction 
coordinate crosses the barrier region. When the electronic "uncertainty" time of charge fluc- 
tuations is shorter than the nuclear dynamics, the transition is adiabatic with the overall 
dynamics evolving on an adiabatic ground-state surface. For small splitting between elec- 



tronic states A « 0.1 — 1, this ground-state adiabatic sur 



ace is often characterized by a 



cusp-shaped potential, j]} [1 U, Q| In fact, it is often argued ?| that majority of natural ET 
reactions are in what we term weakly adiabatic regime, when the reaction is still adiabatic 
but the barrier is quite sharp and characterized by the barrier frequency u a roughly by an 
order of magnitude higher than the medium relaxation frequency ujq. 
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REACTION COORDINATE 



FIG. 1: Schematic energy profiles of reactant (R) and product (P) states of the electron transfer 
reaction coordinate. The reorganization energy E r is the sum of the reaction energy AE and the 
optical excitation energy E* . A stands for energy separation between the energy surfaces due to 
electronic coupling of R and P states. Thermal ET occurs at nuclear configurations characteristic 
to the intersection of the parabolas. 

The dynamics of the reaction coordinate in either of the potential wells can be estimated 
by use of a generalized Langevin equation with a friction term mimicking dielectric response 



of the medium. We are thus 



eft with a standard model of a Brownian particle in a (generally 
In the case of a cusp-shaped potential, a particle approaching 



most surely be pulled towards the other 
by the reciprocal mean first passage 



time-dependent) medium. 0, ^ 
the top of a barrier with positive velocity will a 
minimum. The kinetic rate is then determined \ 
time (MFPT) to cross the barrier wall with positive velocity x > and in a leading order 
in the barrier height AE/ksT yields the standard transition state theory (TST) result 

k cus P = exp (_ A £/^ T ). (!) 
2n 

As discussed by Talkner and Braun, the result holds also for non-Markovian processes 
with memory friction satisfying the fluctuation-dissipation theorem. The TST formula fol- 
lows from the Kramers rate for the spatial-diffusion-controlled rate of escape at moderate 
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and strong frictions r\ 

Here a; a stays for the positive- valued angular frequency of the unstable state at the barrier, 
and o>o is an angular frequency of the metastable state at x — R. For strong friction, rj >> u a 
the above formula leads to a common Kramers result 

k R ^ P = ^±exp(-AE/k B T), (3) 
Znr] 

and reproduces the TST result eq. $1} after letting the barrier frequency tend to infinity 
u a — > oo with ?7 held fixed. Moreover, as pointed out by Northrup and Hynes, [a] in a weakly 
adiabatic case, the barrier "point" is a negligible fraction of the high energy barrier region, 
so that the rate can be influenced by medium relaxation in the wells. The full rate constant 
for a symmetric reaction is then postulated in the form 

k WA = (i + 2k a WA /k D y 1 k a WA (4) 

where k^r A « k cusp and ko is the well (solvent or medium) relaxation rate constant which 
for a harmonic potential and a high barrier (AE > 5k B T) within 10% of accuracy simplifies 
to 

k B T \irk B T J 

In the above equation, the diffusion constant D is related [3J] via linear response theory to 
the longitudinal dielectric relaxation time tl and for symmetric ET reaction {rn^ujQ = 2E r ) 
reads 



k B T k B T 
m LJoT L 2E r r L 

Existence of a well defined rate constant for a chemical reaction requires that the re- 
laxation time of all the degrees of freedom involved in the transformation, other than the 
reaction coordinate, is fast relative to motion along the reaction coordinate. If the separa- 
tion of time scales were not present, the rate coefficient would have a significant frequency 
dependence reflecting various modes of relaxation. Such a situation can be expected in 
complex media, like non-Debeye solvents or proteins, where there are many different types 
of degrees of freedom with different scales of relaxation. Although in these cases the rate 
"constant" can no longer be defined, the overall electron transfer can be described in terms 
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of the mean escape time that takes into account noisy character of a potential surface along 
with thermal fluctuations. 

The time effect of the surroundings ("environmental noises"), expressed by different time 
constants for polarization and the dielectric reorganization has not been so far explored in 
detail, except in photochemical reaction centers a] where molecular dynamics studies have 
shown that the slow components of the energy gap fluctuations are, most likely, responsible 
for the observed nonexponential kinetics of the primary ET process. The latter will be 
assumed here to influence the activation free energy of the reaction AE = E r /4 and will be 
envisioned as a barrier alternating processes. In consequence, even small variations 5E in 
AE can greatly modulate the escape kinetics (passage from reactants' to products' state) 
in the system. If the barrier fluctuates extremely slowly, the mean first passage time to the 
top of the barrier is dominated by those realizations for which the barrier starts in a higher 
position and thus becomes very long. The barrier is then essentially quasistatic throughout 
the process. At the other extreme, in the case of rapidly fluctuating barrier, the mean first 
passage time is determined by the "average barrier". For some particular correlation time 
of the barrier fluctuations, it can happen however, that the mean kinetic rate of the process 
exhibits an extremumjlll. I12I ["K^ that is a signature of resonant tuning of the system in 
response to the external noise. 

In this contribution we present analytical and numerical results for the escape kinetics 
over a fluctuating cusp for the high temperature model ET system with a harmonic potential 
subject to dichotomous fluctuations. We perform our investigations of the average escape 
time as a function of the correlation rate of the dichotomous noise in the barrier height at 
fixed temperatures. In particular, we examine variability of the mean first passage time 



for different types of barrier switching 



ff 



14 Il5| when the barrier changes either between 



the "on-off" position, flips between a barrier or a well, or it varies between two different 
heights. By use of a Monte Carlo procedure we determine probability density function (pdf) 
of escape time in the system and investigate the degree of nonexponential behavior in the 
decay of a primary state located in the reactants' well. 
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II. GENERIC MODEL SYSTEM 

At high temperatures it is permissible to treat the low frequency medium modes classi- 
cally. The medium coordinates are continuous and it is useful to draw a one-dimensional 
schematic representation of the system (Fig. ^) with the reaction proceeding almost exclu- 
sively at the intersection energy. As a model of the reaction coordinate kinetics, we have 
considered an overdamped Brownian particle moving in a potential field between absorbing 
and reflecting boundaries in the presence of noise that modulates the barrier height. The 
evolution of the reaction coordinate x(t) is described in terms of the Langevin equation 



Here £(£) is a Gaussian process with zero mean and correlation < £(£)£(s) >= S(t — s) (i.e. 
the Gaussian white noise arising from the heat bath of temperature T), rj(t) stands for a 
dichotomous (not necessarily symmetric) noise taking on one of two possible values a± and 
prime means differentiation over x. Correlation time of the dichotomous process has been 
set up to ^ with 7 expressing the flipping frequency of the barrier fluctuations. Both noises 
are assumed to be statistically independent, i.e. < £(t)r)(s) >= 0. Equivalent to eq. (J7J) is a 
set of the Fokker-Planck equations describing evolution of the probability density of finding 
the particle at time t at the position x subject to the force —V±(x) = —V'(x) + a±g(x) 



In the above equations time has dimension of [length] 2 /energy due to a friction con- 
stant that has been "absorbed" in the time variable. We are assuming absorbing boundary 
condition at x = L and a reflecting boundary at x = 





d t P(x,a ± ,t) = d x [V£(x)+Td x ] P(x,a±,t) 
- 7P(x, a±, t) + 7P(x, a T , t). 



(8) 



P(L,a ± ,t) = 0, 



(9) 



[V±(x)+Td x ] P{x,a±,t)\ 



x=0 — 



0. 



(10) 
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The initial condition 



P(x,a+,0) = P(x,a_,0) = ~6(x) (11) 

expresses equal choice to start with any of the two configurations of the barrier. The quantity 
of interest is the mean first passage time 



oo L 

MFPT = J dt J [P(x,a+,t) + P(x,a-,t)]dx 



= r + (0) + r_(0) (12) 

with r + and r_ being MFPT for (+) and (— ) configurations, respectively. MFPTs r + and 
r_ fulfill the set of backward Kolmogorov equations 13( 

1 dV±{x) dr±{x) d?r±{x) 
~2 = " 7T±(X) + 1T ^ X) Tx dx~ + dx^ (13) 

with the boundary conditions (cf. eq. (jOJ and (fTTH) ) 



ri(x)| x . =0 = 0, T± (x)\ x=L = 0. (14) 

Although the solution of (fT3j) is usually unique, [3] a closed, "ready to use" analytical 
formula for the MFPT can be obtained only for the simplest cases of the potentials (piece- 
wise linear). More complex cases, like even piecewise parabolic potential V± result in an 
intricate form of the solution to eq. Other situations require either use of approxima- 

n n hn 

tion schemes. |17j perturbative approach \12\ or direct numerical evaluation methods. [181 Il9| 
In order to examine MFPT for various potentials a modified program)^ applying general 
shooting methods has been used. Part of the mathematical software has been obtained from 
the Netlib library. 
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III. SOLUTION AND RESULTS 



Equivalent to equation (fdj) is a set of equations 



du(x) 
dx 
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where new variables have been introduced 
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Since u does not enter the right-hand side of any of the above equations, the syste 
further converted to 
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i.e. it has a form of 



df(x) 

vf 



A(x)f(x)+(3(x). 



A unique solution to |T8|) exists [la] and reads 



f(x) = exp { J A(x')dx' } f(0) 



- exp { j A(x')dx' > J exp ^ - J A(x")dx" } (3{x')dx' 
o Jo I 

A(x)/(0) + A(x)B(x) = A(x)/(0) + C(x) 



with boundary conditions leading to 



/(0) 



r + - r_ 





f(L) 





p(L) 
q{L) 



S 



MFPT is the quantity of interest 



r = r (0) = u(0), 



which can be obtained from 



p(x) = 



du(x) 
dx 



and 



J p(x)dx = u(L) - u(0) = - «(0) = -«(0), 



with 



For the parabolic potential = — l^-(x) = -^r = 2E r x 2 the above procedu 



to 



1 

u u L2T 

2j 2Hx n 

T L 2 T U 



A(x)dx = X 



1 
% 

T L 2 T U 



and 



L L 

t = J J $(x,y)dxdy, 




where 



L 2 H x[<p(x) - 2] LH[<p(L) - 2] 4 7 L 4 T + H 2 x 2 cp(x 



p(x) 



A[H 2 + -fTL 2 ip(L)] 



p(x) 



+ 



+ 



+ 



^pJTm)H x^ 



x H^L 



> y[<p(y)-t\ 

A\m + 1TI?<p(L)\ Jpjxj\ ' " p(y) 
H 2 L^y[<p(y) - 2] x[<p(x) - 2] H 2 y^y) x£(x) 



My) P(x) ATy/pJyjy/pJxj 

4-fL 4 T + H 2 y 2 ip(y) 4 7 L 4 T + H 2 x 2 <p{x) 



ATp{y) 



p(x) 



x 9(y-x), 



p(x) = H 2 x 2 + 2 7 L 4 T, 
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(p(x) = exp 



£(x) = exp 
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FIG. 2: MFPT as a function of the correlation rate of the dichotomous noise for parabolic potential 
barriers switching between different heights H±. Full lines: analytical results; symbols stay for the 
results from MC simulations with At = 10 -5 and ensemble of N = 10 4 trajectories. For simplicity, 
the parametrization L = T = 1 has been used. 

Fig. El displays calculated MFPT as a function of switching frequency 7. Analytical 
solutions are presented along with results from Monte Carlo simulations of eq. (JJJ). As 
a choice for various configurations of the potential, we have probed H± = ±8T; H + = 
8T, H = and H + = 8T, H_ = 4T that set up reorganization energies (2E r = H) and 
heights of barrier in the problem of interest. 

The distinctive characteristics of resonant activation is observed with the average escape 
time initially decreasing, reaching a minimum value, and then increasing as a function of 
switching frequency 7. At slow dynamics of the barrier height, i.e. for values of the rate 7 
less than t+\ _ the average escape time approaches the value r = (r_ — t + )/2 predicted_ by 



theory I 



17 



and observed in experimental investigations of resonant activation. 
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For fast dynamics of the barrier height the average escape time reaches the value associated 
with an effective potential characterized by an average barrier. In comparison to the "on-off" 
switching of the barrier, the region of resonant activation flattens for dichotomic flipping 
between the barrier and a well, and in the case of the Bier-Astumian model (H + = 8T, H = 
AT) when the barrier changes its height between two different values. The resonant frequency 
shifts from the lowest value for the Bier-Astumian model to higher values for the "on-off" 
and the "barrier-well" scenarios, respectively. This observation is in agreement with former 
studies^] aimed to discriminate between characteristic features of resonant activation for 
models with "up-down" configurations of the barrier and models with the "up" configuration 
but fluctuating between different heights. The "up-down" switching of the barrier heights 
produces shorter MFPT and in consequence, higher value of crossing rates for resonant 
frequencies than two other models of barrier switching. 

For each of the above situations we have evaluated probability density function for first 
escape times. Pdfs have been obtained as a result of MC simulations on N = 10 4 tra- 
jectories by use of histograms and kernel density estimation methods. In the resonant 
activation regime, pdf of escape times has an exponential slope, that suggests that the 
reactants' population follows preferentially the kinetics through the state with the lowest 
barrier. Similarly, the exponential decay times of reactants' are observed in the high fre- 
quency limit (7 ~ 10 9 ), when the system experiences an effective potential with an average 
barrier. 11, U, lh, 12, 12, 22, 3,Q 

Apparent nonexponential decay of the initial population is observed at low frequencies 
(7 « 1CT 6 ), when the flipping rate becomes less than tZ 1 . As it is clearly demonstrated in 
Fig. Eland summarized in Table 1, passage times at low frequencies are roughly characterized 
by two distinct time scales that correspond to r_ (different for various switching barrier 
scenarios - cf. first three rows of Table 1) and r + (bottom three rows of Table 1), respectively. 
The inset of Fig El shows pdf zoomed for part of the main plot, where the differences between 
various, model dependent r_ values can be well distinguished. 



IV. SUMMARY 

In the foregoing sections we have considered the thermally activated process that can describe 
classical ET kinetics. The regime of dynamical disorder where fluctuations of the environ- 
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FIG. 3: Semilog plot of the pdf for first escape times in the system. The differences between 
lines fitted to the slopes reflect various r_ values for a given type of barrier switching: (+) for 
H+ = 8T, H_ = -8T; (x) for H + = 8T, H_ = and (*) for H+ = 8T, iT_ = 4T (c/. Table 1). 



H+ H- fitted value static barrier value 

8T -8T 0.09 ±0.01 0.12 

t_ 8T 0.46 ±0.01 0.50 

8T 4T 4.51 ± 0.05 3.43 

8T -8T 68.73 ±0.89 63.01 

t + 8T 71.33 ±3.04 63.01 

8T 4T 79.11 ± 1.98 63.01 

TABLE I: Relaxation time for parabolic potential barriers switching between different heights H± 
at low frequencies. 

ment can interplay with the time scale of the reaction itself is not so well understood. The 
examples where such a physical situation can happen are common to nonequilibrium chem- 
istry and, in particular to ET reactions [^] in photosynthesis. As a toy model for describing 
the nonexponential ET kinetics |5| we have chosen a generic system displaying the resonant 
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activation phenomenon. The reaction coordinate has been coupled to an external noise 
source that can describe polarization and depolarization processes responsible for the height 
of barrier between the reactants and products wells. We have assumed that the driving 
forces for the ET process interconverse at a rate 7 reflecting dynamic changes in the tran- 
sition state. The best tuning of the system and its highest ET rate can be achieved within 
the resonant frequency region. On the other hand, nonexponential ET kinetics can be at- 
tributed to long, time-persisting correlations in barrier configuration that effectively change 
a Poissonian character of escape events demonstrating, in general, multiscale time-decay of 
initial population. 
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